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We analyze the exclusive charmonium J/ip + rjc pair production in e^e~ annihilation using the 
nonfactorized perturbative QCD and the light-front quark model(LFQM) that goes beyond the 
peaking approximation. We effectively include all orders of higher twist terms in the leading order 
of QCD coupling constant and compare our nonfactorized analysis with the usual factorized anal- 
ysis in the calculation of the cross section. We also calculate the quark distribution amplitudes, 
the Gegenbauer moments, and the decay constants for J/ip and rjc mesons using our LFQM. Our 
nonfactorized result enhances the NRQCD result by a factor of 3 ~ 4 at y/s = 10.6 GeV. 



I. INTRODUCTION 

It has been known that the exclusive pair production of 
heavy meson can be reliably predicted within the frame- 
work of perturbative quantum chromodynamics(PQCD), 
since the wave function is well constrained by the nonrela- 
tivistic consideration |l|. However, the large discrepancy 
between the theoretical predictions 0, d, 0, Q based on 
the nonrelativistic QCD (NRQCD) Q factorization ap- 
proach and the experimental results 0, [1] for the ex- 
clusive J/ip + rjc production in e+e~ annihilation at the 
energy -^s = 10.6 GeV has triggered the need of better 
understanding both in the available calculational tools 
and the appreciable relativistic effects. 

A particularly convenient and intuitive framework in 
applying PQCD to exclusive processes is based upon 
the light-front (LF) Fock-state decomposition of hadronic 
state. If the PQCD factorization theorem is applicable, 
then the invariant amplitude for exclusive processes fac- 
torizcs into the convolution of the valence quark distri- 
bution amplitudes (DAs) (/)(a;, q^) with the hard scatter- 
ing amplitude Tf/, which is dominated by one-gluon ex- 
change diagrams at leading order of QCD coupling con- 
stant as ■ To implement the factorization theorem at high 
momentum transfer, the hadronic wave function plays an 
important role linking between the long distance nonper- 
turbative QCD and the short distance PQCD. In the LF 
framework, the valence quark DA is computed from the 
valence LF wave function ^'„(xi,kxi) of the hadron at 
equal LF time r = t + z/c which is the probability ampli- 
tude to find n constituents(quarks,antiquarks, and glu- 
ons) with LF momenta ki — (xi,k_Li) in a hadron. Here, 
Xi and kj^i are the LF longitudinal momentum fraction 
and the transverse momenta of the «th constituent in the 
n-particle Fock-state, respectively. 

The NRQCD factorization approach [3, i, i, Q for 
charmonium production assumes that the constituents 
are sufBciently nonrelativistic so that the relative mo- 
tion of valence quarks can be neglected inside the me- 
son. In this case, the quark DA becomes the 5 function, 
i.e. 4){x^q^) ~ 5{x — l/2)(the so-called peaking approx- 
imation). However, the cross section value 0, H fl 0| 
estimated within the NRQCD factorization approach in 
the leading order of as underestimates the experimental 



data 0, [1] by an order of magnitude. In order to reduce 
the discrepancy between theory and experiment, the au- 
thors in Refs. [1, [l^, [ll|, [l3| considered a rather broad 
quark DA instead of 5-shaped quark DA. 

However, as pointed out in Refs. [3, [3, if the quark 
DA is not an exact 5 function, i.e. k_L in the soft bound 
state LF wave function can play a significant role, the 
factorization theorem is no longer applicable. To go be- 
yond the peaking approximation, the invariant amplitude 
should be expressed in terms of the LF wave function 
^(xi, kxi) rather than the quark DA. In Refs. [l3,[i3lj we 
discussed the validity issue of peaking approximation for 
the heavy pseudoscalar meson pair production processes 
such as e+e" P + P{P = Bs, B, D, D^) using 
the LF model wave function 'i'{xi, kj^.;) cx exp(— Mq //3^), 
where Mq is the invariant mass of the constituent quark 
and antiquark defined by M§ = J2i(^'± + "i^D/^i ^^'^ 
(3 is the gaussian parameter. The gaussian parameter [3 
in our model wave function was found to be related to 
the transverse momentum via (3 = i/(k2_). This relat ion 
naturally explains the zero-binding energy limit as the 
zero transverse momentum, i.e. (Mq) = (mi -I- 7712)^ and 
Xi = rrii/M for (3 — Q. We also found that the heavy 
quark DA is sensitive to the value of (3 and indeed quite 
different from the 5-type DA according to our LFQM 
based on the variational principle for the QCD-motivated 
Hamiltonian [Tsl . [T6| . In going beyond the peaking ap- 
proximation, we stressed a consistency of the formula- 
tion by keeping the transverse momentum kx both in the 
wave function part and the hard scattering part together 
before doing any integration in the amplitude. Similar 
consideration has also been made in the recent investi- 
gation of the relativistic and bound state effects [TtI not 
based on the light-front dynamics(LFD) but including 
the relativistic effects up to the second order of the rel- 
ative quark velocity, i.e. (w^). Such non-factorized anal- 
ysis should be distinguished from the factorized analy- 
sis P, [13, where the transverse momenta are seper- 
ately integrated out in the wave function part and in the 
hard scattering part. Even if the used LF wave functions 
lead to the similar shapes of DAs, it is apparent that the 
predictions for the cross sections of heavy meson pro- 
ductions would be different between the factorized and 
non-factorized analyses. 
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In this work, we extend our previous works [13, ll^ 
pseudoscalar meson pair production to the case of pseu- 
doscalar and vector meson productions and calculate the 
cross section for e+e" J /"ip + r]c process at leading or- 
der of as including effectively all orders of higher twist 
terms. As noted in our results for the quark DA of 
J 1^ and rjc are quite different from the (5-type function. 
We find that the non-factorized form of the form factor 
enhances the cross section of NRQCD result by a factor 
of 3 ~ 4 at ^/s = 10.6 GeV while it reduces that of the 
factorized formulation by 20%. Since the cross section 
for e+e" — > J/ip + rjc is found to be very sensitive to the 
behavior of the end points (a: ^ and 1) in the quark 
DA, we also examine the results of the decay constants 
or equivalently the Gegenbauer moments of J /ip and rjc 
mesons. Since the perturbative corrections of order as to 
the production amplitude has already been obtained 
increasing the cross section significantly, it is important 
to consider the more accurate assessment of cross section 
at the leading order of as ■ 

The paper is organized as follows. In Sec. II, we 
describe the formulation of our light-front quark model 
(LFQM), which has been quite successful in describing 
the static and non-static properties of the pseudoscalar 
and vector mesons [l^ [l^ . The formulae for the quark 
DA, decay constants, Gegenbauer and ^(= xi — X2) nio- 
ments are also given in this section. In Sec. Ill, the 
transverse momentum dependent hard scattering ampli- 
tude and the form factor for e^e" — > J/tp + fJc transition 
are given in leading order of as- The form factors both 
in the factorized and nonfactorized formulations are ex- 
plicitly given in this section. We also show in this section 
that our peaking approximation(i.e. NRQCD) re&Tilt co- 
incides with the one derived from Ma and Si In 
Sec. IV, we present the numerical results for the decay 
constants, quark DAs, Gegenbauer and ^ moments for 
the J/ip and rjc mesons and compare them with other 
theoretical model predictions in addition to the avail- 
able experimental data. The numerical results for the 
e~^e~ — > J/tp + rjc cross section are obtained and com- 
pared with the data [7, 8]. Summary and conclusions 
follow in Sec. V. In the Appendices A and B, we summa- 
rize our results for the helicity contributions to the hard 
scattering amplitudes and the form factor, respectively. 



II. MODEL DESCRIPTION 

In our LFQM [ll,[lll, the momentum space light-front 
wave function of the ground state pseudoscalar and vec- 
tor mesons is given by 



*100 (^i'^^ii'-^i) = '^X^X.,{^i^'^i±)(t>R{Xi,'iii±), (1) 



where (f>R{xi,\<.i±) is the radial wave function and Ti-\'^\^ 
is the spin-orbit wave function obtained by the interac- 
tion independent Melosh transformation from the ordi- 
nary equal-time static spin-orbit wave function assigned 



by the quantum numbers J^'~' . The model wave func- 
tion in Eq. ([1]) is represented by the Lorentz-invariant 
variables, Xi — /P^, k^x = Pi± — XiP± and Xi, where 
and Xi are the momenta and the helicities of con- 
stituent quarks, respectively, and P^^ = {P^ , P^ ,P±) — 
{P° + P^, {AP + Pl_)/P+,P±) is the momentum of the 
meson M. 

The covariant forms of the spin-orbit wave functions 
for pseudoscalar and vector mesons are respectively given 
by 
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where e^{Jz) is the polarization vectors of the vector 
meson, Mq = (k^ + m?)/xiX2 is the invariant meson 

mass square, and X^AiAa '^\i\l'^XiX2 ~ ^ both pseu- 
doscalar and vector mesons. Using the four- vectors pi , p2 
given in terms of the LP relative momentum variables 
(x, kj^) as 



p^^XiP^, P2^X2P^, 

Pi_L = a;iP^ + k_L, P2±=a;2Pj 



(3) 



we obtain the explicit forms of spin-orbit wave func- 
tions for pseudoscalar and vector mesons with the 
longitudinal(e(0)) and transverse(e(-|-l)) polarizations as 
follows 
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where C = ^/2xiX2Mo. For the radial wave function 
(j)R, we use the same Gaussian wave function for both 
pseudoscalar and vector mesons 



4^3/4 
^3/2 



^exp(-fcV2/3^), 



(7) 



where f3 is the variational parameter. When the longi- 
tudinal component k^ is defined by k^ — {x — l/2)A/o, 
the Jacobian of the variable transformation {x,k^} — s- 
k = {k±,kz) is given by dk^/dx = Mo/(4a;iX2). Also, 
the normalization factor in Eq. ([7|) is obtained from the 
total wave function normalization given by 



dx 



f^\M,(J^{x,k^,XrX2)f^l. (8) 
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The quark distribution amplitude (DA) of a hadron in our 
LFQM can be obtained from the hadronic wave function 
by integrating out the transverse momenta of the quarks 
in the hadron, 



167r' 



vI/^3;f(x,ki,AiA2), (9) 



where /i denotes the separation scale between the per- 
turbative and nonperturbative regimes. The dependence 
on the scale is then given by the QCD evolution equa- 
tion 19] and can be calculated perturbatively. However, 
the distribution amplitudes at a certain low scale can be 
obtained by the necessary nonperturbative input from 
LFQM. The presence of the damping Gaussian factor in 
our LFQM allows us to perform the integral up to infin- 
ity without loss of accuracy. The quark DAs for rjc and 
J 1^ mesons are constrained by 



where the decay constant is defined as 



(10) 



mn^qW) = (n) 

for a ?7c meson and 

m-i''q\J/^[P,h)) = fj,^Mj,^e^^{h), 

(12) 

for a J/-0 meson with longitudinal (/i = 0) and 
transverse(/i = ±1) polarizations, respectively. The con- 
straint of Eq. (fTO]) must be independent of cut-off /it up 
to corrections of order A^/u^, where A is some typical 
hadronic scale (< 1 GeV) [1^. For the nonperturbative 
valence wave function given by Eq. ([7]), we take /i ~ rric 
as an optimal scale for our LFQM description of J/?/; and 
r]c- 

The explicit form of the rjc decay constant is given 
by M 



I dx 
2^/& Jo 



167r3 + k^ 



=(/)i^,(a:,ki). (13) 



The decay constants for the longitudinally and trans- 
versely polarized J/i/; meson are given by [20| 
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167r3 + k^ 



respectively. While the constant /, 



Mo + 2m 

(15) 

is known from the 



experiment, the constant /jy^ is not that easily accessible 
in experiment but can be estimated theoretically. 



We may also redefine the quark DA as $^^(j/^)(a;) 
(2\/6//rj^(j/^))0(a;) for the normalization given by 



'^Vc{J/'ip)i^)dx = 1. 



(16) 



The quark DA <i>(x) evolved in the leading order of Q;s(/i) 
is usually expanded in Gegenbauer polynomials Ci?^^ as 



$(a;,/i) = $as(a;) 



l + J2an{^^)Cf/'{2x~l] 



:(17) 



where $as(a;) = 6x{l — x) is the asymptotic DA and the 
coefficients a„(/x) are Gegenbauer moments 19]. The 
Gegenbauer moments with n > describe how much 
the DAs deviate from the asymptotic one. In addition 
to the Gegenbauer moments, we can also define the ex- 
pectation value of the longitudinal momentum, so-called 
^-moments: 

(D = J\cm ^ dxC^ix), (18) 

where $(x) = 2l>(2a; — 1) normalized by (^°) = 1. 

The ^ moments are related to the Gegenbauer moments 
as follows (up to n = 6): 



(C^) = TTT + "2 ;;;; + a4 TT^TTT + "6 
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III. HARD CONTRIBUTIONS TO 

e+e" -» J/iP + ric PROCESS 

For the exclusive process 

e+e- -> 7*(g) ^ J/^iPv) + VciPp), (20) 

the form factor is defined as 

{J/^P{Pv,h)7^,{Pp)\JU0) = e^^-'P-etPvpPp.Tiq^), 

(21) 

where e*{Pv,h) is the polarization vector of the vector 
meson with four momentum Py and helicity h. The cross 
section can be calculated as 

^(e+e- ^ J/V.7?,) = (1 - (22) 

where we neglect the small mass difference between J/tp 
and rjc, i.e. Mh ~ Mj/^ « M^^. 

At leading order of as, the contribution to the form 
factor comes from four Feynman diagrams; one of them is 
shown in Fig.[T] To obtain the timelike form factor T{q^) 
for the process e+e~ — > 7* — > J/^r]c, we first calculate 
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PQCD contribution in the framework of LFD. The sec- 
ondary quark-antiquark pair production can occur only 
through the gluon momentum transfer due to the ra- 
tional light-front energy-momentum dispersion relation. 
The quark-antiquark pair production from the vacuum 
is suppressed in the LFD and the absence of zero-mode 
contribution can be shown by the direct power counting 
method that we presented in our previous work of weak 
transition form factors between pseudoscalar and vector 
mesons [2l|. Taking the perturbative kernel of the BS 
equation as a part of hard scattering amplitude Th, one 
thus obtains 



FIG. 1: One of the four Feynman diagrams for the amplitude. 



the radiative decay process ?7c(-P) + J/'4'{P') 
using the Drell-Yan-West((7+ = + g'^ = 0) frame 



P 



Mi, . + qi 



(23) 



where the four momentum transfer is spacelike, i.e. = 
q^q~ — Qj^ = — < 0- We then analytically continue 
the spacelike form factor !F{q^) to the timelike > 
region by changing — q5_ to q^ in the form factor. 

In the calculations of the form factor !F{q'^), we use the 
'-I- '-component of currents and the transverse(/i = ±1) 
polarization for J/^p given by 



^(±1) 



71 (V.^ 



2q^ 



(24) 



where q^ = qx — i%- For the longitudinal (/i — 0) po- 
larization, it is hard to extract the form factor since 
both sides of Eq. (PT|) vanish for any value. In the 
energy region where PQCD is applicable, the hadronic 
matrix element {J / fp\Jcm\'nc) can be calculated within 
the leading order PQCD by means of a homogeneous 
Bethe-Salpeter(BS) equation for the meson wave func- 
tion. Formally, one may consider a contribution even at 
lower order without any gluon exchange as often called 
Feynman mechanism. However, we don't need to take 
this contribution into account in this work because we 
are considering the production process of heavy mesons 
that ought to require high momentum transfer between 
the primary quark-antiquark pair production and the sec- 
ondary quark-antiquark pair production in order to get 
the final state heavy mesons. Since the final bound state 
wavefunctions satisfy the BS type iterative bound state 
equation, one gluon exchange can be generated by it- 
eration from the wavefunction part even if the scatter- 
ing amplitude formally has no gluon exchange. We thus 
generate the hard gluon exchange from the iteration of 
bound state wavefunction and consider the leading order 



(J/^|J+nhc) - 



E 

A, A' 

xTnix, k^; y, U; qj_; A, A') 
xvI/™o(x,k^,A'), (25) 



where [(Pk] = dx(P]i±/16Tr^ and Th contains all two- 
particle irreducible amplitudes for ^* +qq qq from the 
iteration of the LFQM wave function with the BS kernel. 
On the other hand, the right-hand-side of Eq. (|2ip for 
the matrix element of J"*" is obtained as 



{jmjLivc) 



(26) 



Here, we set = 1 without any loss of generality. 
Therefore, we get the form factor as 



where we combined the spin-orbit wave function into the 
original Th to form a new Ty, i.e. 



V2 



Th = ^y2n\]{,iy,l^)TH{x,k^;y,lr,cir,\,X') 



A, A 
>00 



x7^",^,^(x,k^). 



(28) 



We should note that the form factor J-{q'^) has a di- 
mension [1/GeV] so that the cross section has a dimen- 
sion of [barn] where 1 GeV~^=0.39 mb in the natural 
unit(?i = c = 1). Since the measure [d^k] has the dimen- 
sion of [GeV^] and our radial wave function (f>ij has the 
dimension of [1/GeV], the amplitude Th has the dimen- 
sion of [l/GeV^]. 

The leading order light-front time-ordered diagrams for 
the meson form factor are shown in Fig. [21 where the 
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FIG. 2: Leading order(in as) light-front time-ordered diagrams of the hard scattering amplitude for ??c(-P) — * 7* (9) + J/^{P') 
process. 



energy denominators are given by 
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M 
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2 (U - 2;2q±)^ + (l_L-y2qj 



yi y2 

Dg - D^, Dw = Ds,Dn = Dj, = D2. (29) 

According to the rules of light-front perturbation the- 
ory, the hard scattering amplitude for the diagram Ai in 
Fig. [2] is given by 



Ta, = (-1) 



9{k+) 

z + n^n u{ki+q)-f+u{ki) 
kIDiD2 



4nasCFu{lih^u{ki + q)d%^^v{k2)-l''v{l2) 

(30) 



f(fc„ ) 



ktDiD2 



where Na^ = u{h)-f^'u{kl + q)d%'\{k2)i''v{l2) is the 
gluon-fermion vertex part with the gauge dependent po- 
larization sum d^v for the gluon, Cp = 4/3 is the color 
factor and the notation u{p) denotes actually u{p)/ ^/p^ 
for the internal fermions in a scattering amplitude. In 
the Feynman gauge, the polarization sum d^'^ equals to 
g'^". In the light-front gauge rj ■ A = j4+ = 0, 
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A=l,2 



e^{kg,X)e„{kg,X) 



k^ 



(31) 



where kg ■ e — rj ■ e = 0. Since the gluon propagator 
has an instantaneous part(r/^r/'^/(fc+)^ in the light-front 
gauge), we absorb this instantaneous contribution into 
the regular propagator by replacing kg hy kg = (2/2 ^ 
X2, fcg ,2/2q± +kj_ - lj_), where k' = + - - k^ 
includes the instantaneous contribution. 
Similarly, the hard scattering amplitude for the diagram 
A2 in Fig. O is given by 



(-2)(4^a,CF) 



ktDsDi 



■Na,, 



(32) 



where Na2 has the same form as Nai but with different 
gluon momentum kg = {x2 - y2, kg' - y2qj_ - 
The diagram A2 has also the instantaneous contribution 
and absorb it into k~ = P~ + q~ — (fci -I- q)~ — I2 ■ 
The hard scattering amplitude for the diagram A^ in 
Fig. O is given by 



Ta, = (-2)(4^a,CF) 



k+D^De 



Na„ 



(33) 



where Nas has also the same form as Na^ or Na2- How- 
ever, in this case we have just regular gluon propagator 
with the four momentum kg — {x2 — 2/2, ^g^i Ij- ~ 2/2q± ^ 
kj^) and k~ — k^^//c+ since the diagram does not have 
the instantaneous part. 

Likewise, one can obtain the hard scattering amplitudes 
corresponding to the diagrams Bi. 



6 



If one includes the higher twist effects such as intrinsic 
transverse momenta and the quark masses, the LF gauge 
part proportional to leads to a singularity although 
the Feynman gauge part gives the regular amplitude. 
This is due to the gauge-invariant structure of the ampli- 
tudes. The covariant derivative I?^ = dfj, + igA^ makes 
both the intrinsic transverse momenta, kj^ and and 
the transverse gauge degree of freedom gA± be of the 
same order, indicating the need of the higher Fock state 
contributions to ensure the gauge invariance (2^ . How- 
ever, we can show that the sum of six diagrams for the 
LF gauge part(l/fc+ terms) vanishes in the limit that the 
LF energy differences A^, and Aj, go to zero, where A^, 
and Ay are given by 

XiX2 

yiy2 " 

Details of the proof can be found in our previous 
work [l^ . In this work, we follow the same procedure 
presented in Ref. and calculate the higher twist ef- 
fects in the limit of A^, = A^ = to avoid the involve- 
ment of the higher Fock state contributions. Our limit 
Aj; = Aj, = 0(but (k\) = P ^ Q) may be considered 
as a zeroth order approximation in the expansion of a 
scattering amplitude. That is, the scattering amplitude 
Th may be expanded in terms of LF energy difference 
AasTH = [Tff](o) + A[rff](i) + A2[r^](2) where 
[Th-](°) corresponds to the amplitude in the zeroth order 
of A. This approximation should be distinguished from 



the zero-binding(or peaking) approximation that corre- 
sponds to Af = mi + 1712 and k_L = f3 = 0. The point of 
this distinction is to note that [Tf/]^°) includes the bind- 
ing energy effect(i.e. kj^,lj^ ^ 0) that was neglected in 
the peaking approximation. 

In zeroth order of A, one can show that the energy 
denominators entering in the diagrams Ai and Bi(i — 
1,2,3) have the following relations 

V2 + V4 = Vi,V2+V5^0,V2 + Vs^V'r, (35) 

where Vi = limA^=A„=o£'i and 

= -(a;2qi-F2ki-qi)/a;i, 
= -(z/2qi-2U-qi)/yi, 
2?2 = --[xjyjqi + xjll + ylkl ~ 2x2y2{x2l± ■ q± 
-y2^± • q_L + k_L • U) + (y2 - X2 fm'^] 

/[X2y2{y2 - X2)]. (36) 

As one can see from Eqs. ^ and (O, the leading 
twist(LT)(i.e. neglecting tranverse momenta k^ and 1_l) 
helicity contributions for ?7c(/')+ 7* (q) — > J/V'(P', e(+l)) 
process come from two leading helicity AH =^ — 
A,,^| = 1 components, i.e. ti~*TT and it^TT as in the 
nonrelativistic spin case, i.e. -^(| ti) ^ I it)) for ijc 
meson to | IT) for J/ip meson. 

In Table HJ we summarize our results for the hard 
scattering amplitude Th for these two leading helicity 
AH = 1 components in zeroth order of binding energy 

limit. For instance, J^Li ^a/'"^^^ and J^Li ^s/^^^^ 
are obtained as 
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(37) 



To derive the final results in Eq. p7p , we use the following for the diagram A and 

identities obtained from Eqs. (I29p and ([35| : ^ ^ ^ ^ 
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- v^v^v,. - v,v,(-K, ,39, 

~ I?i I?2 ' '"^^'^ diagram B. The above identities lead to 0{y2 — 
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TABLE I: Leading helicity contributions to the hard scattering amplitudes, y^Ti^TT) T^^^^^^'. The momentum variable 
pRW represents p^^^^ = ± ipy 



\AH = 1\ 


Na^ = 


Na2 = 




= Ay 


= 0) 




Nb2 = Nbs 






= Ay 


= 0) 


n^TT 

iT^TT 










-TT) 


2m{xi I 


-yik-xiy2q)^ 








-TT) 


2 in ( X 2 


X2yiy2 

l-y2k-x2y2q)^ 


{y2-X2)'Dl'D2 


-TT) 


2m{x2 


a;2I/l!/2 

l-y2k-x2y2q)^ 


(y2 

i 




-^TT) 




x-,yTy9. 


(y7-x-?.)T>-,T>2 






XTy-,y2 


(1/2 


-^2)D2Dr 





2:^2) + 6'(a;2 - ya) = 1 in Eq. 

By adding all six LF time-ordered diagrams, we obtain 



(n^TT)i(o) 



3 

E 

i=l 



T 



(Ti^TT) 



T 



(Ti^TT) 
Bi 



(40) 



A=0 



Similarly, one can easily obtain the helicity it^TT contri- 
bution ti 
Table H 



bution to the hard scattering amplitude, yj^J^T^TT) fj.^^ 



In the numerical calculations for the higher twist con- 
tributions, one may keep effectively only the leading or- 
der of higher twist terms such as k^/q5_, l^/q^, and 
k± due to the fact that <C q^^ and 1^ <C 

in lar ge rn omentum transfer region where PQCD is appli- 
cable [13, H^- As shown in our previous work J4J, this 
can be done by neglecting the subleading higher twist 
terms accordingly both in the energy denominators and 
the numerators for the hard scattering amplitude Th- 
This procedure is very similar to the recent investigation 
of the relativistic and bound state effects not based on 
the LFD but including the relativistic effects up to the 
second order of the relative quark velocity, i.e. 
Indeed, our numerical result neglecting the higher or- 
ders of k^/q5_, Ij^/qi and • l±/q^ is very close to 
that presented in Rcf. 17] (see Section IV). However, in 
this work, we include all higher orders of k^ /q^ , l^^/q^ 
and k^ • l±/q\. This corresponds to keep effectively all 
higher orders of the relative quark velocity beyond (w^). 
We compare our full result with the one neglecting the 
corrections of order 0{{v'^)). 



Using Eq. (|28p . we then obtain the leading helicity 
contributions to the hard scattering amplitude combined 
with the spin-orbit wave function in zeroth order of A as 
follows 



(0) 



V2^iit 
qL '^TT 



(Ti^TT)TpOO 



H 



(iT^TT)-r>00 



H 









m 


m - 


May+2m 



a/to2 + k^ y^m^TYf 



M 



H, 



(41) 



where 
Mh 



KSTTasCprn 



1 1 



xiX2yiy2{y2 - X2)q'^ 
-{xiyi + X2y2)k^q" 



+ 



xiX2yiy2{y2 - X2) 
y2{xl+xl) 



V1V2 V2VS 
X2{xiyi + X2y2) 



V1V2 



(42) 



and l^q^ = 1± • q± + x qj^|. Accordingly, the lead- 
ing helicity contributions to the form factor lead to the 
following non-factorized form 



Qc I [d^k][dH]Mx,'^±WHf UR{yA 



+Qc{mi ^ TO2) 
2QcJ [d^k][d^l]- 
1 



m 



2^ki 
1^ 



bR{x,'ki_)MH 



Oy 



2m 



0i?(i/,U), 



(43) 



where Qc = 2/3 is the charge fraction of charm quark in 
the unit of e. 

In the leading twist (LT) limit neglecting the transverse 
momenta, the hard scattering amplitude in Eq. (I42|) is 
reduced to 

= 2 —[xi{xiyi+X2y2) + yi{xl+xl)]. 44 

XiX'^ywiq^ 

Then, the form factor in Eq. ([^5]) factorizcs into the 
convolution of the nonpcrturbative valence quark DAs 
(t>rj^(j/^){x, ^) with the perturbative hard scattering am- 
plitude M^J^: 



2Qc dx dy(j)y^^ {x , ^x) 
Jo Jo 

xM^{x„y^,q'^)<l)^j/^{y,n), (45) 



where (lr^^^{y,^) is the quark DA for the transversely 
polarized J/ip meson. Furthermore, in nonrelativistic 
QCD(NRQCD) limit(i.e. peaking approximation) where 
the longitudinal momentum fractions are given by Xi = 



8 



?/i « 1/2 with M/i « 2mc, the quark DAs for both 77c and 
J /ip mesons become (5-type functions, i.e. 



^(a;-l/2), (46) 



^,.iJMy-^-' 2V6 

and the form factor in NRQCD Hmit is reduced to 



„^ /r,e/j/V' 2567ras(M)CF 



(2^6)2 
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(47) 



where the superscript 5 for the quark DA in Eq. 
and the form factor in Eq. (1471) represents the NRQCD 
resuh. Our NRQCD result [4^ is exactly the same as 
that derived from Ma and Si in Ref. [3 (see Eqs. (16) 
and (21) in 

Other subleading helicity contributions to the hard 
scattering amplitude that show up as next-to-leading or- 
der in transverse momenta are summarized in Tables IIVI 
and fVl of the Appendix A. 



IV. NUMERICAL RESULTS 

In our numerical calculations, we use our LFQM [Tsl . 
[l6j parameters {mdicc) obtained from the meson spec- 
troscopy with the variational principle for the QCD mo- 
tivated effective Hamiltonian. In our LFQM, we have 
used the two interaction potentials Vqq for rjc and J/'0 
mesons: (1) Coulomb plus harmonic oscillator (HO) po- 
tential, and (2) Coulomb plus linear confining potential. 
In addition, the hyperfine interaction essential for the 
distinction between J/tjj and rjc mesons is included for 
both cases (1) and (2), viz.. 



conf-|^ + ^^S^Q-V^(r), (48) 



where Vconf = ^f"^ for the HO potential and br for the lin- 
ear confining potential, respectively. For the linear con- 
fining potential, we use the string tension 5 = 0.18 GeV^, 
which is rather well known from other quark-model anal- 
ysis commensurate with the Regge phenomenology 24]. 
The other potential model parameters arc then fixed by 
the variational principle for the central Hamiltonian with 
respect to the Gaussian parameter (3. For instance, the 
model parameters for the linear confining potential are 
obtained as a = —0.724 GeV, rUc = 1-8 GeV, and the 
strong couphng constant as(M) = 0.313 defined by 



as(/i) 



127r 



(33-27V/)ln(^VA2)' 



(49) 



where Nf = 4 is the number of active flavors(w, d, s and 
c). At scale ji ~ rric for charmonium, our value of A = 
162 MeV, the scale asscociated with nonperturbative ef- 
fects involving light quarks and gluons, is consistent with 
the usual Kqcd — 200 MeV. Our value of as = 0.313 is 
also quite comparable with other quark model predictions 



TABLE II: Decay constants [MeV] of rjc and J/4> obtained 
from our variational parameters (mc = 1.8, /3) [GeV] and com- 
pared with the experimental data. 





Linear 


HO 


HO' 


Exp. 




(/3 = 0.6509) {P 


= 0.6998) {P 


= 0.7278) 




fva 


326 


354 


370 


335 ± 75 [29] 


fj/i> 


360 


395 


416 


416 ± 6 [30] 


fj/i, 


343 


375 


393 





such as 0.35, 0.45, 0.30 ~ 0.38, and 0.314 from ISGW2 
model 2^i, Corn ell p otential model [2^, Bodwin-Kang- 
Lee(BKL) model [23] (in next-to- leading order in Ug), and 
relativistic quark model [iTj, respectively. Lattice mea- 
surements of the heavy-quark potential yield the values 
for effective coupling as of 0.22 in the quenched case and 
approximately 0.26 in the unquenched case (28j . The HO 
potential model parameters are obtained in a similar way 
as in the case of the linear potential. We should note 
that the root-mean-square value of the transverse mo- 
mentum in our LFQM is equal to the Gaussian (3 value, 
i.e. ^/l^JTc = Pec- 

In Table [Tll we summarize our results for decay con- 
stants of rjc and J /^p obtained from our variational pa- 
rameters (mc = 1.8 GeV, Pec =0.6509 GeV) for the 
linear potential [second column] and (toc = 1-8 GeV, 
Pec =0.6998 GeV) for the HO potential[third column]. 
Our results for the decay constants fn^ = 326 [354] 
MeV, fji^ = 360 [395] MeV, and /J^^ = 343 [375] MeV 
obtained from the linear [HO] potential parameters are 
quite comparable with the current experimental data, 
(/r^Jexo = 335 ± 75 MeV ^ and (/j/^)exp = 416 ± 6 
MeV [20] as well as other theoretical model calculations 
such as the QCD sum rules [3l|, [13] in which the decay 



constants were obtained as 



346 MeV, /. 



412 



MeV and /J/,/^ = 409 MeV. As a sensitivity check of our 



variational parameters, we include in Table [TTl another 
Gaussian parameters (toc = 1.8 GeV, Pec =0.7278 GeV) 
to fit the central value of the experimental J /ijj decay 
constant. We denote this as HO' in Table HIl In the fol- 
lowing numerical calulations, we present all of these three 
cases (linear, HO, HO') to show the parameter sensitivity 
of our results. 

The shape of the quark DA which depends on (rric. Pec) 
values is important to the calculation of the cross section 
for the heavy meson pair production in e^e~ annihila- 
tions. We thus show in Fig. [3] the normalized quark DA 
for Tje and ipccix) — (t>nai^) ~ 4'j/^{x) obtained from 
linear (dotted line), HO (solid line), and HO' (dashed line) 
potentials compared with the ones obtained from Bon- 
dar and Chernyak(BC) ^ (dot-dashed line) and from 
QCD sum rules 31](doublcdot-dashed line). As one can 
see from Fig. [31 our quark DA 4>ecix) obtained at scale 
fi ~ rUe practically vanishes in the regions x < 0.1 and 
X > 0.9 where the motion of cc pair is expected to be 
highly relativistic. However, our results for quark DA are 
certaintly wider than the delta function-type(i.e. Pec 
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FIG. 3: The leading twist distribution amplitudes (^(a;) for r]c 
and J/tp[(j)ric{x) ~ (j}.J/ii{^)] obtained from our LFQM com- 
pared with the ones from Bondar-Chernyak model and 
QCD sum rules [H. 



limit) NRQCD results 0, 13, H, which do not take into 
account the relative motion of valence quark-antiquark 
pair. Our results also show that the shape of quark 
DA becomes broader and more enhanced at the endpoint 
region(a; — > or 1) as the Gaussian parameter /5(or equiv- 
alently transverse kj^-size) increases. In comparison with 
other theoretical model calculations, we find that our re- 
sult is quite consistent with the one obtained from QCD 
sum rules [sij at scale ~ rric but much narrower than 
the one obtained from BC @. As will be discussed later, 
the cross section for double-charm production is indeed 
very sensitive to the end point behavior of the quark DA. 

In Table IIIIl we list the calculated (^") moments up 
to n = 6 for the rjc and J/4' DAs at scale /i ~ tUc and 
compare with other model estimates. Our central, up- 
per and lower values are obtained from HO, HO' and 
linear parameters, respectively. Since the ^ moments 
for J/ip meson with the longitudinal polarization are al- 
most the same as those with the transverse polarization, 
our results imply that = = (C")t, which 

is also confirmed by the recent QCD sum rule calcula- 
tions [31I. [s^l. Furthermore, the ^ moments between r^c 
and J/ijj mesons are not much different from each other as 
one can see from Table iHll Our results for the ^ moments 
are in good agreement with those obtained from other po- 
tential models [H, as well as QCD sum-rules [3l|, [33, 
but disagree with the predictions obtained from BC [91 
and BKL ^ models. While NRQCD predictions [34 1 
for the second and fourth moments are in agreement with 
our model but disagree for the higher moment (^^). This 



disagreement for the moment (^^) may be ascribed to the 
end point behavior(i.e. relativistic correction) of quark 
DA. 

Within our model calculation, the relative quark ve- 
locity can be obtained from the relation ~ 

2'\J + ~ 2mc in the center of mass frame, where 
rrir is the reduced mass. From this relation, we obtain 



777,2 ^ 2/3V?772,i.e. 



= 0.30 



-0.02 

-0.04' 



(50) 



where the central, upper, and lower values are from the 
HO, HO', and linear potential parameters, respectively. 
Using the dimensional regularization at leading order of 
Q!s, the authors in Refs. [3l|, [s^l derived the relation be- 
tween the relative velocity of quark-antiquark pair in- 
side the charmonium and the ^ moment as {C^)n^ ~ 
(e),7/^ = (7;")/(77-f 1)(77 = 2,4,6). Applying this for- 
mula to our model calculations, we get the relative as 
{v^)cc = 0.251q'q2i where again the central, upper, and 
lower values are from the HO, HO', and linear potential 
parameters, respectively. The results obtained from our 
LFQM and QCD sum rule methods are not only in an 
agreement with each other but also quite consistent with 
the value {v^)cc « 0.3 used in NRQCD [3,[23. Note that 
one gets the quark DA (j){x) ^ 5{x — 1/2) in the Hmit 
w — > while (j){x) ~ (l>as{x) = Qx{l — x) as v ^ 1. As 
noted in QCD sum rule calculations [3l|, the moments 
(^") are proportional to according to the NRQCD v- 
scaling rules {rricV^)^ <C {rricv)'^ <C tti^. It is not diffi- 
cult to see that the ^-moments obtained from our LFQM 
and QCD sum rules [3l| satisfy these rules. However, as 
discussed in [3l|, the BC moments ^9,] break the NRQCD 
w-scaling rules and the quark DA obtained from [§| cor- 
responds to the QCD sum rule result [3l[ defined at scale 
/i ~ 10 GeV rather than at /x ~ rric- In our LFQM cal- 
culation, we would overestimate the experimental values 
of decay constants for J/ip and rjc if we were to use the 
shape of BC distribution to get the cross section value 
of double-charm production consistent with the experi- 
mental data. Thus, it seems misleading to claim that the 
cross section of e+e^ J/ip + 77c in [9| led to a good 
agreement with the experiment. 

The corresponding Gegenbauer moments obtained 
from Eq. (jl9p are given by 



a2{^l ^ nic) 

04 (/I ~ 777c) 

ae{fi ~ 777c) 
for T^c meson and 

a2(M ~ 777c) 
ai{fj. ~ 777c) 

ae{fj. ~ 777c) 



-0 S39+° °" 

0.00»_g Q21, 

0-082:i:o;o22' 
+0 0027+° °'^''^ 

-0.0112' 



-0 343+0 011 

U.O^O_Q Q20, 

0.087:f:o;Q2o, 
-0 001 ^s+o oo^s 



(51) 



(52) 



for J/tp meson, respectively. Since a„ for J/ip with longi- 
tudinal polarization and with transverse polarization 
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TABLE III: The S, moments {C")>7c s^nd {C")j/i/> — ~ for r/c and J/ip distribution amplitudes obtained from our 

LFQM at the scale fi ~ rric and compared with other model({^"),,^ « estimates. Our central, upper, and lower values 

are obtained from the HO, HO', and linear potential parameters, respectively. 



(D 


Ours 


Ours 


BuchmuUer Cornell 


BC [9] 


BKL [27] 


NRQCD [34] 


QCD sum 








Tye model [33] model [26] 








rules [31, 32] 


n = 2 


0.084±ii;rr 


0.082iHli^ 


0.086 0.084 


0.13 


0.019 


0.075 ±0.011 


0.070 ± 0.007 


n = 4 


'^■"-'J^ ' -0.003 


coieit;;!!;!^ 


0.020 0.019 


0.040 


0.0083 


0.010 ±0.003 


0.012 ±0.002 


n = 6 




o.oo46±;j:;ii!^^ 


0.0066 0.0066 


0.018 


0.0026 


0.0017 ± 0.0007 0.0031 ± 0.0008 



are not much different from each other, we do not distin- 
guish them in our model calculation. 

In Fig. m we show s^JF(s) for e+e^ — * J/V' + ?7c pro- 
cess. The left panel of Fig. |4] shows the results obtained 
from the central value /3 = 0.6998 GeV of our model 
parameters displaying different (leading and subleading) 
helicity contributions. The right panel of Fig. 0] shows 
the sensitivity of our model predictions with all helic- 
ity contributions when the gaussian model parameter (3 
changes as shown in Fig. [31 In the left panel, the dot- 
ted and short-dashed lines represent the results obtained 
from the non-relativistic peaking approximation J-^{s)[ 
Eq. (|47p ] and the leading twist(LT) factorized form fac- 
tor T^"^ {s)[ Eq. HSl)] taking into account the relative 
motion of valence quarks, respectively. The long-dashed 
line represents the higher twist (HT) nonfactorized form 
factor J^^'^(s)[Eq. (|^ ] obtained by including the trans- 
verse momenta(kj^, Ij^) both in the wave function and 
the hard scattering part. Note that JF* (s) (dotted line), 
.F'""'"(s) (short-dashed line), and ^^^"""(5) (long-dashed line) 
are obtained from the leading helicity contributions. The 
solid line represents our full solution •?^(Ait=o+AH=i)('^) 
including all(leading plus subleading) helicity contribu- 
tions summarized in the Appendix A (Tables IIVI and W\ . 
Among the subleading helicity contributions, we find 
that only tT^TT(dot-dashed line) and (Ti^Ti) + (iT^iT 
) (double-dot-dashed line) helicity contributions give a 
sizeable effects and other subleading helicity contribu- 
tions are negligible. As shown in Fig.[31 we find that while 
•^^•^(AH=o+AH=i)('*) about 2 times larger than s'^T^{s) 
but 10% smaller than s'^T^^{s) at ^/i = 10.6 GeV. It 
is also interesting to note that our s^-^(Ait=o+AH=i) 
takes over s^T^'^[s) for ^/s > 13 GeV region, although 
s^T^'^{s) with leading helicity components approaches 
to s'^!F^^{s) as s ^ 00. As one can see from Fig. |4l 



form factor 

1 ^-2 



obtained from our calculation shows 
as s — » 00 which is the expected QCD 



the 

scaling behavior [l^, [H, [H, [13, IH, [s^ for the transi- 
tion form factor between pseudoscalar (0 '") and vector 
(1 ) mesons. 

In Fig. 5, we show leading order in as contribution to 
the cross section for e+e" J /ip + i]c- The left panel of 
Fig. [5] shows the results with leading and subleading he- 
licity contributions using the HO model parameters. The 
right panel of Fig. [5] shows the sensitivity of our model 
predictions with all helicity contributions when the gaus- 
sian model parameter /3 changes as shown in Figs. [3] and 



m The line codes are the same as in Fig. U) As one can 
see from the left panel of Fig. 5, our peaking approxima- 
tion result (dotted line) is consistent with the previous 
NRQCD estimates in Refs. which is an order of 

magnitude smaller than the experimental data 0, Q • We 
should note from the left panel of of Fig. 5 that our higher 
twist result (solid line) including all helicity contributions 
enhances the peaking approximation result by a factor of 
3 ~ 4 at ^/s = 10.6 GeV while it reduces that of the lead- 
ing twist result by 20%. As discussed in Section III, our 
higher twist results (ctht) include all orders of k^/q^, 
ll/qi and • l±/q^ to keep effectively all higher or- 
ders of the relative quark velocity beyond {v'^). If we keep 
only the leading order of these terms (k^/q5_, lj_/qj_ and 
k_L • l±/q^), our results would correspond to include the 
relativistic effects up to the order of ('y^)[l3]- Our pre- 
dictions for the cross section at a/s = 10.6 GeV obtained 
from peaking approximation((T5), leading twist(CTLT) and 
higher twist (ctht) are given by 



= 2.34 



-0.50 
-0.69 

10.57 



[fb]. 



= 8.76 



+ 1.61 



-2.84 



[fb], (53) 



where the central, upper and lower values are obtained 
from HO, HO' and linear potential parameters, respec- 
tively. Our prediction of a'^^!^~'^^'^^~^^ {J/ 4' + Vc) 
reduces by about 10% from the value in Ea. ((53| to 
7.68l;2'g6[ft'] when we keep only the leading order of 
k^/q'^, 1^/q^ and • l±/q^. It is interesting to note 
that our reduced value 7.68I2 eel^] indeed very close 
to the result 7.8 [fb] obtained in the recent investigation 
including the relativistic effects up to (u^)[l3|- 

As a sensitivity check, we show in Fig. [5] the parame- 
ter (mc, f3) dependence of the cross section for e~^e~ 
J/ 4' + rjc using the nonfactorized higher twist form factor 
with all helicty contributions. We also show in Fig. [H] 
the decay constants corresponding to the end point mass 
values, nic = 1.4 GeV and 1.8 GeV. The cross section 
increases as f3{mc) increases (decreases). As one can also 
see from Fig. [HI the cross section is more sensitive to the 
variation of the gaussian parameter than to the variation 
of the charm quark mass. 

The experimental results are 

aiJ/ip + T]c) X B"^" [> 2] = (25.6 ± 2.8 ± 3.4) [fb] , (54) 
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FIG. 4: The form factor s^T{s) for e^e~ — *■ J/Tp+Vc- The dotted, short-dashed, long-dashed and soUd Une in left panel represent 
the peaking, leading twist (LT), higher twist (HT) results with the leading helicty contributions and HT one with all helicity 
contributions, respectively. The dot-dashed and double-dot-dashed line represent the dominant subleading contributions. The 
right panel represents the HT results including all helicity contributions obtained from linear(dotted line), HO (solid line), and 
HO'(dashed line) model parameters, respectively. 



by Belle (filled circle in FigO and 

(t( J/V + Vc) X S"^ [> 2] = (17.6 ± 2.8t ^-j) [fb] , (55) 

by Babar [8] (filled square in FigO, where i?''<=[> 2] is the 
branching fraction for rjc decay into at least two charged 
particles. Considering an enhancement by the factor of 
1.8 from the corrections of next-to- leading order (NLO) 
of as , it might be conceivable to raise our leading as 

order result cr'^^~^^^^~^\j /ip + rjc) in Eg. ([55)1 by this 
factor and get a value close to the above Babar data. 
However, it would be necessary to make detailed NLO 
investigation within the LF PQCD framework before we 
can make any firm conclusion. 



V. SUMMARY AND CONCLUSION 

We investigated the transverse momentum effect on 
the exclusive charmonium J/ip + rjc pair production in 
e^e~ annihilation using the nonfactorized PQCD and 
LFQM that goes beyond the peaking approximation. 

Our LFQM calculation based on the variational prin- 
ciple for the QCD-motivated Hamiltonian ^15, 16] shows 
that the quark DAs for J /ip and rjc take substantially 
broad shape which is quite different from the i5-type DA. 
If the quark DA is not an exact 5 function, i.e. the 
relative motion of valence quarks can play a significant 
role, the factorization theorem is no longer applicable. 



In going beyond the peaking approximation, we stressed 
a consistency by keeping the transverse momentum kj^ 
both in the wave function part and the hard scattering 
part simultaneously before doing any integration in the 
amplitude. Such non-factorized analysis should be dis- 
tinguished from the factorized analysis where the trans- 
verse momenta are seperately integrated out in the wave 
function part and in the hard scattering part. Even if 
the used LF wave functions lead to the similar shapes of 
DAs, predictions for the cross sections of double-charm 
productions are apparently different between the factor- 
ized and non-factorized analyses. We found that the 
higher twist contributions including all helicity contri- 
butions enhanced NRQCD result by a factor of 3 4 at 
= 10.6 GeV while it reduced that of the leading twist 
result by 20%. We also found that the cross section for 
e+e" J/^'Vc process at y/s = 10.6 GeV is more sen- 
sitive to the variation of the gaussian parameter than to 
that of the charm quark mass. Our results showed that 
the cross section increases as /3(toc) increases(decreases). 

In conclusion, LFQM/PQCD analysis showed that the 
relativistic correction(i.e. non-delta function) of the 
light-front wave function is very important to under- 
stand the large discrepancy between the NRQCD re- 
sult and the experimental data given by Eqs. ([5^ and 
([55)1 . While there have been considerations of broaden- 
ing the quark DA to reduce the discrepancy between the 
theory at the leading order of Us and the experimental 
results [1, [13, [iH, a recent calculation of corrections of 
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FIG. 5: The cross section for e^e J /ij^rjc with leading and subleading helicity contributions using the HO model parame- 
ters(left panel) and with all helicity contributions using the linear, HO, and HO' model parameters(right panel). 
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FIG. 6: The parameter (mc,/3) dependence of the cross sec- 
tion for e~ J /ijj + rjc- 



next-to-leading order (NLO) of as leads to an enhance- 



ment of the theoretical prediction by the factor about 
1.8 [3]. This factor may enhance our result in the lead- 
ing order of as to fit the current experimental results. 
However, more detailed investigation is necessary prior 
to any firm conclusion on this issue. 
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APPENDIX A: HELICITY CONTRIBUTIONS TO 
THE HARD SCATTERING AMPLITUDE 

In this appcnix A, we summarize the helicity contribu- 
tions (Ai, A2) (^^1,^2) to the hard scattering amplitude 



H 



for the r]c{P) — > 'y*{q) + J/ip{P') process. 
In Tables HVl and fVl we summarize our results for the 
helicity contributions to the hard scattering amplitudes 
Ta and Tb for the diagrams in Fig. [21 where 
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TABLE IV: Helicity contributions to the hard scattering amphtude Ta in Fig. (2] 



lAi^l Hchcitics 

TT^TT 

J, -j, -1- -1' 



Na, 



Na, 



Na, 



J:^TaA.A. = Ay ^ 0) 



' (1/2-3:2) ^ I'iX'2 

' (V2-X2) ^ •Clf 2 

(a2-3:2)E'lE'2 

(M-a:2)E'ir'2 
STrasCj. , FA+iG 

' (1/2-2:2) ^ I'l-E'2^ 
' (1/2-2:2) ^ P1P2 







Pai - iC 



Pa2 + iQA 
Pa2 - iQA 



Pa3 + iQA 
Pa3 - iQA 



(y2- 



2m {yi-xi)(y2-X2) 

xi^2yiy2 

2m'^ (yi-xi)(y2-X2) 

^\^2y\y2 
Fai + iGA 

Fai - iGA 



N 



(n^iT) 

Ai 



r(iT^Ti) 

Fa2 + iGA 
Fa2 — iGA 



(n^iT) 

Ai 

(iT^Ti) 
Ai 

Faz + iGA 
Fa3 — iGA 



2 1 



■X2)V2^ 

■n) 



(y2- 
(Ti- 
Ai 

(iT^n) 
1 1 

2 1 

(i/2-3:2)'P2 ' 



(H2- 



n^TT 

it^TT 

TT^n 

TT^iT 



2m(^xil — yik — X2yiq) 

^2yiy2 
2m(x2l — y2k — X2y2q)' 

xiyiy2 

2m(x2l — y2k-X2y2q)'" 

a:iSlS2 
2m(xil — yik — X2yiq)^ 

a:2!/l!/2 
2m(xil — yik — X2yiq)'^ 

xiX2y2 
2m{x2l — y2k-X2y2q)^ 
xiX2yi 
2m{x2l — y2k — X2y2q)'^ 

xix2yi 
2m(xi l~yik-X2yiq)^ 
xiX2y2 



N 



Ai 



Ai 
■{IT 
Ai 
■{IT 



Ai 



Ai 



Ai 



Ai 



TT) 




TT) 










{y2-X2)T>iT>2 


ID 




ii) 










{y2-X2)'Dl'D2 


TT) 




TT) 


{y2-^2)'Dl'D2 








ii) 


N^A^ 


ii) 


STrag C p 








(y2-^2)T>lT>2 


Ti) 




Ti) 


STTd^Cp 








(y2-X2)'DiT'2 


Ti) 




Ti) 


SttcisCp 


7\7"(J-J- 






{y2-X2)'Di'D2 


TT) 




iT) 


SwctsCp 


< 






(y2-X2)ViT>2 


XT) 


7V7-(-l-l— * 

'-^Ai 


iT) 


SnasCp 


AT" (ii 
^1 






(y2-^2)T>lT>2 



♦TT) 
♦ii) 
-TT) 
♦ii) 
♦Ti) 
♦Ti) 
♦iT) 
♦iT) 



2 


TT^ii 
















ii^TT 















TABLE V: Helicity contributions to the hard scattering amplitude Tb in Fig. [2] 



|A_ff| Helicities 

TT^TT 



Ne 



N, 



B2 



S3 



0) 







Pbi + iQs 
Pbi — iQs 

2m^(yi-xi){y2-X2) 

xiX2yiy2 
2m (yi-xi)(y2-X2) 

xix2yiy2 

Fbi + iGs 
Fbi - iGb 



Pb2 + 
Pb2 - 



N 



(Ti- 
Sl 



iQs 
iQB 

HT) 



r{iT 
'Si 

Fb2 + 
FS2 - 



^Ti) 



iGb 
iGB 



Pb3 + iQs 

Pbs - iQB 

Ar(Ti-iT) 
^^Si 

A7(iT^Ti) 

^^Si 

Fb'A + iGB 
Fb3 — iGB 



(H2 



y2-X2) I- 'D2T>g 
irasCp i Pb-^Qb 
X2) I- 15208 

(m-^2)'D2'o» 

Sttocs Cp 
iy2-^2)'D2'Ds 

(1/2-3:2) ^ 'E'21'8 

(y2-X2) I- 'P2P8 



(j/2 



-^2)T>5' 
2 1 



(a2 
(Ti- 



"'Si 

^(iT- 



a:2)I'5J 
iT) 



'Si 



Ti) 
2 



(!/2 



- 3:2)05 J 
2 1 



(l/2-^2)P5 



Ti^TT 
it^TT 



2m(3:ii — 1/1 fc -3:1^29) 

a:2!/l!/2 
2m(3:2i — 1/2^-3:21/29)^ 
a:iI/lI/2 
_ 2m(3:2i — 1/2^-3:2^29) 

!i:iai!/2 
_ 2m(3:ii -1/1 fc -3:11/2 !?)-"• 

!i:2!/l!/2 
_ 2m(3:ii — i/ifc — 3:ii<2g) 

a:ia;2!/2 
_ 2m(3:2t — 1/2^-3:2 2/2 9)"'" 
xiX2yi 
2m(x2l — y2k — X2y2q)'^ 

xi^2yi 
2m(xxl — yik-xiy2q)^ 
xiX2y2 



r(Ti^ 


TT) 




TT) 


8-ITOlsCp 


^^Si 


Si 






(y2-^2)'D2'Ds 


r{Ti^ 


ii) 




ii) 


SlTCsCp 


^^Si 


Si 






(y2-^2)'D2'Ds 


r{iT-^ 


TT) 


AtUT- 

J^Si 


TT) 


Sttocs Cp 


< 


Si 






(V2-^2)'D2'Ds 


r{iT^ 


ii) 




ii) 


SttCsCp 


< 


Si 






(y2-^2)'D2'Ds 


r{TT- 


Ti) 




Ti) 


SttCsCp 




Si 






(V2-^2)'D2'Ds 


r{ii^ 


Ti) 


Ar(ii— * 
^^Si 


Ti) 


SttcisCp 


7vr(ii 
^^Si 


Si 






{y2-^2)'D2'Ds 


r(TT^ 


iT) 


^^Si 


iT) 


SttCsCp 


^^Si 


Si 






{y2-^2)'D2'Ds 


r(ii-^ 


iT) 


7vr(ii ' 

^^Si 


iT) 


SttCsCp 


^^Si 


Si 






iy2-^2)'D2'D8 



♦TT) 
♦ii) 
♦TT) 
♦ii) 
♦Ti) 
♦Ti) 
♦iT) 
♦iT) 



2 
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Pai = 

Qa = 

Fai = 

Fa2 = 

Ga = 



-2 



x\viy2C^i_ + xiX2\\ + yi?/2ki + X2{xiyi + X2y2)(l± • q±) 



xiX2yiy2 

+2x2yiy2(k± • q±) + {xiyi + X2j/2)(m^ + kj_ • 1±) 



(a^i - yi) 



[D2 + £"4] = 



[xi - yi) 



[D2 + Di\, 



Pa 



{xi - yi) 



[D2 + Di], Pa3^Pa 



{xi - yi) 



[D2~D4], 



xiX2yiy2 

-2 
xiX2yiy2 

Fa 



X2{xi - y2)\i± X q^l - 2xiyiy2|k_L x q_L| + (a;22;2 + xiyi)\'k±_ x 

2 



2 



(xi - yi) 



[D2 + D^], Fa3=Fa 



{xi - yi) 



(xi - yi) 
[D2-Di], 



[D2 + Di] = Fa 



(xi - yi) 



[D2+Di], 



xiX2yiy2 
for the diagrams Ai and 
-2 



2xiyiy2\'k± x q_L| - xsjU x q^| + (1 - 2a;iyi)|k_L x 



(Al) 



Pb2 

Qb 
Fbi 



xiX2y2<i_L + 2:1x21^ + yiy2^± - y2{xiyi + X2y2){'^i^ ■ q±) 



a:ia;2yiy2 

-2xiX2y2(U ■ q_L) + {xiyi + X2y2){rn^ + k± • U) 



= Pr 



{x2 - y2) 



[D^ + D^], Pb3 = Pe 



iy2 - X2) 



(2/2 - X2) 



[Dq + D7] = Pb 



(2/2 - X2) 



xiX2yiy2 
-2 



a;ia;2yiy2 
Fb2 — Fi 



-{xiyi - X2y2){\l± X k^l - j/2|q± X k^l), 

kj_ • 1_L - y2(k_L • q±) + {xiyi + X2y2)m^ 

2 



{x2 - 2/2) 



[Do + £'71, F, 



B3 — £b 



(y2 - X2) 



(2/2 - 2:2) 



[Z?9 + D7] = £s 



{y2 - X2) 



-(2/2|q± X k_L| - |U X k_L|), 



Gb — 

xiX2yiy2 

for the diagrams S^, respectively. 



(A2) 



As an illustration, we show how to obtain the hard as well as the total amplitude Th — Ta + Is for the 
scattering amplitudes Ta ~ Yl^i=i (sixth column in (Ti^Ti) contribution. Using the identities Eqs. 
Table |IV| and Tb = ELi ^is. (sixth column in Table |V| and ^ in Sec.III, we obtain 



En! 



(n-'Ti) 



A=o 



d{y2-x2){PAi+iQA) , e{x2-y2){PA2+iQA) ^0{x2-y2){PA3 + iQA)' 



(2/2 - X2)'Di'D2 



{X2 - y2)V3'Di 



ix2 - y2)V5Ve 



-8na,CF{PA + iQA) 
-IGnasGp 



0{y2 - X2) 9{x2 - 2/2) / 



(2/2 - X2)'Di'D2 {x2 - 2/2) K'D^'Di V^Vq 
9{y2-X2)V2+Vi 0{x2-y2) ( V2+Vi V2-Vi 

{y2-x2Y V{D2 {x2-y2Y\ -D^-Di 



5-^6 



'?>TTasCF{.PA + iQA) 



1 



(2/2 - X2)'Di'D2 



— STTasCp 



(2/2 - X2)^V2 



(A3) 
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and 



E 



(Ti^Tl) 



Bi 



-SnasCp 



A=o 



0{y2-X2){PB^+iQB) , 0{y2-X2){PB3 + iQB) , e{x2 ~ V2){Pb2 + iQ b)' 



'8T:a,CF{PB +iQB) 



(2/2 - X2)'DjI)s 

d{y2 - X2) 



(2/2 - X2)'Dii'Di2 {X2 - y2)'Ds)'Dio 

1 1 \ 6{x2-y2) 



iy2 - X2) \D-jV^ X'iiX'12 / {x2 - y2)'D9'Dia 

e{y2-X2) fVg+Vj ^ Vg^VA , 9{x2 - y2) Vg + Vj 



(2/2 - X2)^ 

SnasCpiPB +iQB) 



1 



(2/2 - X2)V2V, 



VIIV12 



{x2 ~ 2/2)^ 'DgT^io 
2 

iy2 - X2Y'Dz 



(A4) 



where the first terms in Eqs. (jA3j and (|A4ll proportional 
to l/(j/2 — 2:2) and the second terms proportional to 
1/(2/2 — 2^2)^ are related with the Feynman gauge and 
the LF gauge parts, respectively. By adding all six LF 
time-ordered diagrams, we obtain 



E 



T 



(Ti^Ti) 



T 



(Ti^Ti) 

Bi 



(2/2 - X2) 



{Pa + iQ) 



V1V2 



A=0 

{Pb 



'D2T>s 



(A5) 



(1) AH = contributions: 



Tpllt-TpOO 



my/2 

CxCy 



TplltTpOO 



2/1 M( 



y2MQy + m 



M, 



V2 

C'xCy 

V2 



Qy 

m 4 



CxCy 



2m 



■7?llt-7?00 



iL _ TplltTpOO 



Mqj, + 2m 



1 



Moj; + 2m 
(2) AH = 1 contributions: 



(BI) 



i.e. the singular LF gauge parts cancel each other and 
only finite Feynman gauge parts contribute to the am- 
plitude. Similarly, we obtain other helicity contributions 
to the hard scattering amplitude as shown in Tables IIVI 
andlVl 



APPENDIX B: HARD SCATTERING 
AMPLITUDE COMBINED WITH RELATIVISTIC 
SPIN-ORBIT WAVE FUNCTION 



-7?llt-7?00 



Tpllt-TpOO 



V2 



Tpllt-TpOO 



^llt-r>00 _ 



CxCy 



1 



L\2 



V2 

CxCy 

J2_ 

CxCy 



Moy + 2m 
yiMoy + m 



Moy + 2m 
yiMpy + m 



V2 

CxCy 



May + 2m 
2/2 + m 



2m 



V2_ 

CxCy 



2/2M, 



M, 



Oy 



2m 



(B2) 



In this appendix B, we list the leading and subleading 
helicity contributions to the hard scattering amplitude 
combined with the relativistic spin-orbit wave function, 
where the subleading helicity contributions show up as 
next-to-leading order in transverse momenta. That is, 
the subleading helicity contributions vanish at leading 
twist. 

We first consider the relativistic spin-orbit wave func- 
tions for pseudoscalar and vector(with transverse polar- 
ization e = +1) mesons given by Eqs. ^ and re- 
spectively. Besides the leading helicity (in tranverse mo- 
menta) contributions coming from two AH — 1 contri- 
butions(i.e. ti~^TT and it^TT)i the subleading helicity 
contributions are as follows. 



where Cx = \/2xiX2Mox and Cy = ^J2yiy2MQy. Since 
the hard scattering amplitudes vanish for AH — 2 cases, 
we do not consider them here. 

Next, we obtain the hard scattering amplitude com- 
bined with the spin-orbit wave function. 

(1) AH ~ contributions: 



IGTra^Ci^m 



(2/2 - X2)CxCyq'^ [ M^y + 2m 



2/1 



Oy 



PA+iQA , PB+iQB 
V1V2 



V2VS 



(B3) 
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2/2 -^Oj/ + TO 



(y2 - X2)CxCyq'^ 



Moy + 2m 



PA-iQ 



(B4) 



iitT-(TT^TT).^oo 



— 167rasCF 
(2/2 - Xi)CxCyq^ 

'FA + iGA 



TT 



m 



V1V2 



+ 



Moy + 2m 
Fb+iGb 



'D2VS 



(B5) 



(2/2 - X2)CxCyq'^ [Moy + 2m 
Fa-IGa , Fe 



x{k^l^l^q^) 



iGi 



V{D2 



V2VS 



(B6) 



^#(Ti^iT) 



V2 

qL it 



Ti 



327rasCi?m^(j/i — xi) 
xiX2yiy2CxCyq'^ 



y2Moy + m 
Mov + 2m 



x{l^q^) 



+ 



J 00 



327ra5CFm^(j/i — xi) 
xiX2yiy2CxCyq'^ 

x(/V)^ ^ ^ 



yiMpy + m 
Mo„ + 2m 



2?22?8. 



(B8) 



V2 

-327rasC'F?Ti^ 
2:22/12/2(^2 - X2)C^Cyq^ 



1 



Moy + 2m 



xil^q^) 



xi\\ - yil^k^ - X2yil^q^ 



V1V2 



+ 



xil'i - yil^k^ - xiy2l^q^ 
V2VS 



(BIO) 



T#(TT^Ti) 



V^^iitjiCTT^Tij^oo 
327rasCFm 



TT 



XlX2y2{y2 - X2)CxCyq'^ 



yiMoy + m 



Moy + 2m 



x(Z^5«) 



xik^l^ — 2/ik5_ — X2yik^q^ 



V1V2 



+ 



xik^l^ - yik']^ - xiy2k^q^ 



V2V8 



, (Bll) 



T#(ii^n) 



V2 



lltr7i(ii^Ti).^00 



V 'T 



32'jTasCFTn 



xiX2yi{y2 - X2)CxCyq^ 
^^^,^n^\x2l-k^-y2^^- 



yiMoy + m 



Moy + 2m 
- X2y2q 



+ 



V1V2 

X2l^k^ - 2/2ki - X2y2q^k^ 



V2VS 



, (B12) 



(B7) r#(TT->iT) 



V2 

TpiitTi(TT^iT)T700 

■^'nt 'MT 

S2'KasCFm 



xiX2yi(y2 - x2)CxCyq^ 



y2MQy 



M, 



Oy 



2m, 



x(Z^««) 



X2k^l^ - y2k± - X2y2k^q^ 



V1V2 



X2k^l^ - yaki - X2y2k'^q 



L„R 



+ 



V2V8 



, (B13) 



(2) AH = 1 contributions: 



9" 



Ti 



-i2'na,CFm^ 



xiym{y2 -X2)CxCyq'^ 



1 



X{l^q^) 



Moy + 2m 
X2II - y2l^k^ - X2y2l^q^ 



+ 



V1V2 

X2l± - 2/2^^fc^ - X2y2l^q^ 
V2VS 



(B9) 



T#(iWT) 



V2. 



^7^i^t4ii^iT)^oo 
32'jTasCF'm 



XlX2y2{y2 - X2)Ca:Cyq^ 



y2Moy 



m 



Moy + 2m 



x{l^q^) 



xil^k^ - 2/ik5^ - X2yiq^k^ 



+ 



V1V2 

xil^k^ - yik\ - xiy2q^k^ 



2i^8 



• (B14) 
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